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ABSTRACT 

There is observational evidence showing that stellar and solar flares occur in a sim- 
ilar circumstance, although the former are usually much more energetic. It is expected 
that the bombardment by high energy electrons is one of the chief heating processes of 
the flaring atmosphere. In this paper we study how a precipitating electron beam can 
influence the line profiles of Lya, Ha, Can K and A8542. We use a model atmosphere 
of a dMe star and make non-LTE computations taking into account the non-thermal 
collisional rates due to the electron beam. The results show that the four lines can be 
enhanced to different extents. The relative enhancement increases with increasing for- 
mation height of the lines. Varying the energy flux of the electron beam has different 
effects on the four lines. The wings of Lya and Ha become increasingly broad with the 
beam flux; change of the Can K and A8542 lines, however, is most significant in the 
line centre. Varying the electron energy (i.e., the low-energy cut-off for a power law 
beam) has a great influence on the Lya line, but little on the Ha and Can lines. An 
electron beam of higher energy precipitates deeper, thus producing less enhancement 
of the Lya line. The Lya/Ha flux ratio is thus sensitive to the electron energy. 

Key words: line: profiles - stars: activity - stars: atmospheres - stars: flare - stars: 
late-type. 



1 INTRODUCTION 

Stellar flares have been frequently detected in both broad 
band photometry and in spectral lines. A large stellar flare 
can be several orders of magnitude more energetic than its 
solar analogy. During a stellar flare, the continuum emis- 
sion in UV and visible bands is significantly enhanced and 
is usually characterised by a blue colour (e.g. Pagano et al. 
1997); on the other hand, spectral lines show an increased 
net emission, with a magnitude varying from line to line 
(e.g. Catalano & Frasca 1994; Montes et al. 1996) ; lines show 
asymmetries and macroscopic velocity fields of several hun- 
dreds of km s _1 possibly exist (e.g. Houdebine et al. 1993; 
Gunn et al. 1994; Montes et al. 1999), implying that the 
flaring atmosphere is very dynamic. 

A detailed comparison between stellar and solar flares 
was made by Haisch, Strong & Rodono (1991). Despite the 
big difference in magnitude and scale, stellar and solar flares 
still show many similar features in the enhancement of lines 
and continua, implying further that the basic heating pro- 
cesses and emission mechanisms are similar in these two 
phenomena. It is widely accepted that a solar flare occurs 
through the reconnection of magnetic fields, which releases 
a large amount of thermal energy and results in the accel- 



eration of charged particles. Accelerated particles then pre- 
cipitate downwards to heat the lower atmosphere, leading 
to a dynamic evolution of the atmosphere (chromospheric 
evaporation and condensation). Although there is no direct 
information on spatially resolved magnetic activity on stars 
like that observed on the Sun, some observations still show 
evidence that a stellar flare may result from the interaction 
between an old magnetic structure and a new emerging flux 
(Catalano & Frasca 1994). Simon, Linsky & Schiffer (1980) 
proposed a speculative scenario of major long-lived RS CVn 
flares in which the component stars have very large corotat- 
ing flux tubes, which occasionally interact. Recently, Ruben- 
stein & Schaefer (2000) proposed that superflares, which are 
newly detected on F-G main-sequence stars by Schaefer, 
King & Deliyannis (2000) , are caused by magnetic reconnec- 
tion between fields of the primary star and a close-in planet. 
In these cases, the basic flare mechanism is also magnetic 
reconnection. 

Therefore, we can postulate that the eruption of a stel- 
lar flare involves the acceleration of charged particles which 
then bombard the lower atmosphere. This process not only 
provides direct heating through Coulomb collisions, produc- 
ing the bulk mass motions as revealed in the observed spec- 
tra (e.g. Houdebine et al. 1993), but also yields non-thermal 
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collisional excitation and ionization of the ambient atoms. 
The latter increases the source function and thus enhances 
the line and continuum emission. This issue is important in 
the spectroscopy and, in particular, the atmospheric mod- 
elling of stellar flares. It is then desirable to have a quan- 
titative assessment of the non-thermal effects on both the 
line and continuum emission. Numerical computations in- 
cluding the non-thermal effects have been done mostly for 
solar flares (e.g. Fang, Henoux & Gan 1993; Henoux, Fang 
& Gan 1995; Ding & Fang 1997; Ding & Schleicher 1997; 
Ding 1999; Fang, Henoux & Ding 2000). Recently, Ding & 
Fang (2000) studied the role of non-thermal electrons in the 
optical continuum of stellar flares. They found that an elec- 
tron beam with a large energy flux can produce a [/-band 
brightening and a U — B colour that are roughly compara- 
ble with the observed values of a typical large flare. This 
work investigates how the non-thermal electrons influence 
the lines of hydrogen (H I Lya and Ha) and the lines of ion- 
ized calcium (Call K and A8542) which are frequently used 
in the spectroscopy of stellar atmospheres (e.g. Hawley & 
Pettersen 1991; Panagi, Byrne & Houdebine 1991; Houde- 
bine & Doyle 1994a, 1994b; Houdebine & Stempels 1997). 



2 COMPUTATIONAL METHOD 

When an electron beam precipitates downwards from the 
corona, the beam electrons collide elastically with the ambi- 
ent electrons and protons, leading to a direct heating of the 
plasma; meanwhile, the beam electrons collide inelastically 
with the ambient atoms and produce a non-thermal excita- 
tion and ionization of the atoms. As in Ding & Fang (2000) , 
we first evaluate the rate of energy deposition corresponding 
to the latter process, which reads 



-€)n H A'*(*-2)f| ( Wi 



«2 1 d u 

2 + 'V-V 



(1 - u)±+Z 



where £ is the ionization degree of the ambient plasma, K = 
2-7re 4 , A' the Coulomb logarithm for inelastic collisions. We 
assume a power law flux spectrum for the electron beam, 
characterised by a power index, 8, a low-energy cut-off, Ei, 
and a total energy flux, T\ . N is the column density and iVi 
is the depth penetrated by electrons of an energy Ei, which 
is expressed as 
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where is the cosine of the initial pitch angle, taken to 
be unity here. The expression for j3 and 7 can be found in 
Emslie (1978). Moreover, in equation (1), 



1, N>N U 
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(3) 



We adopt an atomic model of four bound levels plus 
continuum for hydrogen and a model of five bound levels 
plus continuum for ionized calcium. The non-thermal exci- 
tation and ionization rates are computed as 



o„$/ni, 



(4) 



where ni is the hydrogen ground level population. The co- 
efficients taken from Fang et al. (1993): in cgs units, 



012 = 2.94 x 10 10 , 013 = 5.35 x 10 9 , aw = 1.91 x 10 9 and 
die = 1.73 x 10 10 for H; an = 2.38 x 10 10 , oi 6 = 4.25 x 10 10 
and ai c = 4.69 x 10 10 for Call. The above formulae for non- 
thermal collisional rates are derived empirically using the 
collisional cross-sections for various transitions (Fang et al. 
1993). It is very convenient to include these rates into the 
computation of model atmospheres. 

We do non-LTE computations using a method similar to 
that described in Gan & Fang (1987) and Fang et al. (1993). 
The general procedure is as follows. Given a model atmo- 
sphere and a prescribed electron beam, we solve iteratively 
the equations of hydrostatic equilibrium, particle conserva- 
tion, radiative transfer and statistical equilibrium, including 
the non-thermal collisional rates (Cy), for both H and Call. 
After the computations reach convergence, we then calculate 
the line source function and the opacity, and finally the line 
profile as 



Fx = 2n 
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Broadening mechanisms include radiative damping, the 
Doppler effect and the Stark effect. Note that for all the 
lines involved, we assume complete frequency redistribution 
(CRD), while the effect of partial frequency redistribution 
(PRD) may be important for some lines (e.g. Lya and Call 
K). However, CRD is still a reasonable assumption in the 
present investigation since our main purpose is to show the 
effect of non-thermal electron beams, not to reproduce ex- 
actly the observed line profiles. 



3 INFLUENCE OF ELECTRON BEAMS ON 
THE LINE EMISSION 

Many authors have investigated how the spectral lines, in- 
cluding the H 1 Lya and Ha, and the Ca 11 K and A8542 lines, 
change with the modifications of atmospheric models of, for 
example, late-type stars (e.g. Linsky et al. 1979; Houde- 
bine & Doyle 1994a, 1994b; Houdebine, Doyle, & Koscielecki 
1995; Houdebine & Stempels 1997; Short & Doyle 1997, 
1998; Mauas 2000). These lines are shown to be good diag- 
nostics of the atmospheric conditions in the late-type stars. 
Generally speaking, the Call K line, in particular, the Kl 
minimum, can be used to determine the temperature struc- 
ture around the temperature minimum region; emission in 
the hydrogen Balmer lines (Ha etc.) is related to the tem- 
perature and temperature gradient in the chromosphere; the 
Lya flux relies sensitively on the position of the transition 
region, i.e., the coronal pressure. However, we will show be- 
low that the presence of a non-thermal electron beam in 
a flare star can also change these spectral features with- 
out modifying the atmospheric model. Therefore, in order 
to construct reliable semi-empirical models or make proper 
spectral diagnostics for these flare stars, we need first to 
study quantitatively the effect of non-thermal electrons on 
the above-mentioned lines. 



3.1 Effect of varying the beam flux 

As in Ding & Fang (2000), we adopt the semi-empirical 
model of the dMe star, AD Leo, proposed by Mauas & 
Falchi (1994), as the base model, and then introduce an 
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Figure 1. Line profiles of Lya, Ha, Can K and A8542 computed 
from a model atmosphere of the dMe star AD Leo (Mauas &: 
Falchi 1994), considering the effect of a precipitating non-thermal 
electron beam. Different line styles show the effect of varying the 
energy flux, T\, of the beam: (solid line), 10 10 (dashed line) 
and 10 12 (dotted line) erg cm - 2 s _1 . fn all cases, the low-energy 
cut-off, Ei, is 20 keV, and the power index, <5, is 3. 



electron beam into the atmosphere to compute the emergent 
line profiles. We fix the model atmosphere by only chang- 
ing the energy flux of the beam, covering a range of 10 9 - 
10 12 erg cm -2 s _1 , which are typical values for a solar flare 
case. Fig. [I] displays the four line profiles in different cases 
of beam fluxes. Generally speaking, the four lines undergo 
a drastic change when the beam flux increases. However, 
there is clearly a difference in the evolutional behaviour of 
these lines. When the beam is not so strong {Ti ^10 
ergcm~ 2 s _1 ), the Lya line shows only an increased emis- 
sion in its far wings, while the inner part, including the emis- 
sion peaks and the line centre, is not increased but some- 
how decreased. When the beam flux further increases, the 
whole line is enhanced significantly, except for the line cen- 
tre, which remains as weak as in the quiescent state. The 
Ha line shows instead a monotonic increase in the whole 
profile with the beam flux. In the case of a strong beam 
(J~i ~ 10 12 erg cm -2 s _1 ), both the Ha and Lya lines be- 
come very broad. Quite differently, change of the Call K 
and A8542 lines is most significant in the line centre. While 
the wings of the K line can also be enhanced to some ex- 
tent, the wings of the A8542 line are slightly reduced by the 
non-thermal effects. 

To show more quantitatively the variation of the lines, 
wo plot in Fig. | the net integrated flux, AF, defined as the 
flux integrated over the whole line profile with the quiescent 
value subtracted, against the electron beam flux. The inte- 
gration extends to AA = ±6 A for Lya and Ha lines and to 
A A = ±2 A for Can K and A8542 lines, respectively. It can 
be seen that the extent of line flux enhancement increases in 
the order of Call A8542, K, Ha and Lya. From line profile 
computations, we can find that this is just the order of line 
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Figure 2. Net integrated fluxes of Lya, Ha, Call K and A8542 
against the energy flux of the electron beam. In all cases, the 
low-energy cut-off, Ei, is 20 keV, and the power index, <5, is 3. 



formation heights, spanning a range from the upper photo- 
sphere or lower chromosphere to the upper chromosphere. 
So the results are compatible with the general picture that 
the electron beam precipitates downwards from the corona. 
We note that the negative value of AF for Lya in weak 
beam cases is due to a reduction of emission in the inner 
part of the line, as described above. 



3.2 Effect of varying the electron energy 

To check the effect of non-thermal electrons of different en- 
ergies, we have also computed the line spectra for electron 
beams with a fixed energy flux but different low-energy cut- 
offs. The line profiles and the net integrated fluxes are plot- 
ted in Figs. |§] and ^ respectively. Fig. q shows that, varying 
Ei does not alter obviously the profiles of Ha, Call K and 
A8542. However, it has a great influence on the Lya profile, 
namely, increasing E\ reduces the non-thermal effects on 
the line intensities. This result is not surprising. As the Lya 
line is formed in the very upper chromosphere, electrons of 
a higher energy penetrate deeper and thus leave less energy 
in upper layers to excite the hydrogen atoms. 

In the above computations, the electron beam is as- 
sumed to originate in the corona, i.e., at the top of the model 
atmosphere. If the beam originates in a lower layer, its effect 
could become more significant for lines formed in the lower 
atmosphere. Such a scenario has been proposed to be a pos- 
sible origin of a flare-like phenomenon on the Sun, Ellerman 
bombs (Ding, Henoux & Fang 1998). 



3.3 Discussion 

To illustrate the different behaviours of the four lines in re- 
sponse to the electron beam, we plot in Fig. |H] the height 
distribution of the line source function in different cases of 
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Figure 3. Line profiles of Lya, Ho, Call K and A8542 computed 
from a model atmosphere of the dMe star AD Leo (Mauas &: 
Falchi 1994), considering the effect of a precipitating non-thermal 
electron beam. Different line styles show the effect of varying the 
low-energy cut-off, E\, of the beam: 20 (dotted line), 60 (dashed 
line) and 100 (solid line) keV. In all cases, the energy flux, T\ , is 
10 12 erg cm - 2 s — 1 , and the power index, <5, is 3. 
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Figure 4. Net integrated fluxes of Lya, Ha, Can K and A8542 
against the low-energy cut-off of the electron beam. In all cases, 
the energy flux, jFi, is 10 12 erg cm - 2 s — 1 , and the power index, 
<5, is 3. 
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Figure 5. Height distribution of the source function of the four 
lines in the atmosphere of AD Leo in the presence of an electron 
beam. The line styles have the same meaning as in Fig. [l]. 



beam fluxes. It shows that the beam can make the source 
function of Lya greatly increased in a broad region from the 
chromosphere to the photosphere. This can be understood 
that the hydrogen ground level is depopulated while the ex- 
cited levels are overpopulated with respect to the beam-free 
case. Accordingly, the Lya wings are significantly enhanced. 
In the very upper layers, however, the source function of Lya 
could be slightly reduced, which results in a decrease of the 
line centre intensity in some cases (see Fig. 0). Similar results 
have been obtained for the Lya line in solar flares (Henoux 
et al. 1995). The source function of Ha is also increased but 
less significantly than Lya. Note that the Stark effect, due 
to an enhanced electron density, plays an important role in 
the broadening of the Ha line. Fig. ^ displays the height 
distribution of the electron density, which is shown to be 
raised by up to three orders of magnitude in the chromo- 
sphere, highlighting the non-thermal ionization effect of the 
beam. 

The increase of the source function of hydrogen lines 
in photospheric layers is mainly caused by a backwarming 
effect, i.e., the enhanced radiation from the chromosphere, 
while the direct penetration of non-thermal electrons has 
less effect. In contrast, the source function of the Can lines 
is mainly increased in the chromosphere; it is still coupled 
to the local Planck function in the photosphere. In addition, 
the Can lines are less affected by the Stark effect. Therefore, 
the enhancement of these lines is most apparent in the line 
centre, as shown in Fig. 

The quantitative results of line and continuum en- 
hancement due to an electron beam are certainly model- 
dependent. In all the above computations, we have assumed 
a cool background atmosphere. In fact, the model atmo- 
sphere should change with the activity level of the star, es- 
pecially when flare occurs. A model of the flaring state of AD 
Leo, for example, was proposed by Mauas & Falchi (1996). 



© 2001 RAS, MNRAS 000,[j]-| 



The role of non-thermal electrons in the hydrogen and calcium lines of stellar flares 5 




10° i i I i i i I i i i I i i i I i i i I i i i I i 
1000 800 600 400 200 
Height (km) 

Figure 6. Height distribution of the electron density in the at- 
mosphere of AD Leo in the presence of an electron beam. The 
line styles have the same meaning as in Fig. [j]. 

Compared with the model in the quiescent state, a flare 
model has a higher temperature in the chromosphere and 
even in the photosphere, and a lower position of the transi- 
tion region (a higher coronal pressure). This model itself can 
produce an enhanced line emission. However, a higher coro- 
nal mass in the flaring atmosphere would consume a part of 
the electron energy, so that the effect of the electron beam 
is less significant than in the case of a quiescent atmosphere. 
The real situation might be a combination of the effects of 
a heated atmosphere and a non-thermal electron beam. 

In addition to the line flux, we find that the Lya/Ha 
flux ratio is in positive dependence on the electron beam 
flux, but, apparently, in negative dependence on the low- 
energy cut-off. In the case of T\ = 10 12 ergcm _2 s _1 , 
Ei — 20 keV and <5 = 3, the ratio amounts to ~ 6. There- 
fore, it seems that an electron beam can produce most spec- 
tral features that are usually ascribed to a modification of 
the model atmosphere. As long as the electron beam exists 
during the flaring process, its effect should be taken into 
account in order to properly construct the semi-empirical 
models, otherwise the chromospheric temperature will be 
overestimated. 



4 CONCLUSIONS 

There is observational evidence showing that a stellar flare 
may occur in a similar circumstance to the eruption of a 
solar flare. A common flare scenario is the reconnection of 
magnetic fields. Although a stellar flare, in particular, the 
flares in late-type stars, can be several orders of magnitude 
more energetic than a solar flare, their emission features, in 
both lines and continua, still share many similarities. It can 
be postulated that the flaring process in a star also involves 
the acceleration of energetic particles (protons or electrons) , 
carrying most energy to heat the lower atmosphere. 

We study in this paper how a precipitating electron 
beam can influence the line profiles of Lya, Ha, Call K 
and A8542. We make non-LTE computations taking into ac- 



count the non-thermal collisional rates due to the electron 
beam. The results show that the four lines can be enhanced 
to different extents. The relative enhancement increases with 
increasing formation height of the lines. Varying the energy 
flux of the electron beam has different effects on the four 
lines. The wings of Lya and Ha become increasingly broad 
with the beam flux; change of the Call K and A8542 lines, 
however, is most significant in the line centre. Varying the 
electron energy (i.e., the low-energy cut-off for a power law 
beam) has a great influence on the Lycv line, but little on 
the Ha and Call lines. An electron beam of higher energy 
precipitates deeper, thus producing less enhancement of the 
Lya line. The Lya/Ha flux ratio is thus sensitive to the elec- 
tron energy. Together with our computations for the contin- 
uum emission (Ding & Fang 2000), these results provide a 
diagnostic tool for the non-thermal processes in stellar flares. 
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